Evolution from BCS Superconductivity to Bose Condensation: Analytic Results for 

the crossover in three dimensions 
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We provide an analytic solution for the mean-field equations and for the relevant physical quan- 
tities at the Gaussian level, in terms of the complete elliptic integrals of the first and second kinds, 
for the crossover problem from BCS superconductivity to Bose-Einstein condensation of a three- 
dimensional system of free fermions interacting via an attractive contact potential at zero tempera- 
ture. This analytic solution enables us to follow the evolution between the two limits in a particularly 
simple and transparent way, as well as to verify the absence of singularities during the evolution. 

PACS numbers: 74.20_z, 74.20.Fg, 67.40.-w 



I. INTRODUCTION 

The interest in the crossover from BCS superconduc- 
tivity to Bose-Einstein (BE) condensation has quite re- 
cently increased both experimentally and theoretically 
after the appearance of new angular resolved photoemis- 
sion data on underdoped cuprate samples, which indicate 
the persistence of a gap in the single-particle excitation 
spectrum at temperatures well above the superconduct- 
ing critical temperature.trH Previous to these data the 
suggestion that superconductivity in cuprates (as well as 
in other "exotic" materials) might require an intermedi- 
ate approach between the two (BCS and BE) limits has 
been emphasized especially by Uemura.Q 

From the theoretical point of view, evolution from BCS 
superconductivity to BE condensation has beep, consid- 
ered originally by Nozieres and Schrnitt-Rink,Ll -follow- 
ing the pioneering works by EaglestJ and Leggett.EI More 
recently, the crossover problem was addressed in Refs. 
|§|-|lo|, prompted by the experimental suggestion that the 
superconducting coherence length in cuprates is consid- 
erably shorter than in conventional superconductors. In 
all these works known results were recovered analytically 
in the two (BCS and BE) limits. However, no analytic 
result was obtained for the most interesting crossover re- 
gion, which had thus to be treated numerically (with the 
sole exception of the two-dimensional case considered in 
Ref. 0). 

In this paper we provide an analytic solution in the 
three dimensional case covering the whole crossover re- 
gion for a number of physical quantities evaluated at the 
mean-field level and with the inclusion of Gaussian fluc- 
tuations, by considering a system of fermions in free space 
at zero temperature mutually interacting via an attrac- 
tive contact potential. For this case we are, in fact, able 
to express these physical quantities in terms of the com- 
plete elliptic integrals of the first and second kinds, whose 
analytic properties and numerical values are extensively 
tabulatcdEj Our solution enables us to interpolate in a 
rigorous fashion between the two (BCS and BE) limits, 
thus avoiding the problems which occur with a full nu- 



merical approach. Although our analytical solution has 
been unavoidably obtained with simplifying assumptions 
(namely, at zero temperature and using a free-particle 
dispersion relation and a contact interaction potential), 
it might nevertheless be regarded as a reference solution 
with which numerical solutions obtained for more com- 
plicated cases could be compared. 

The plan of the paper is as follows. In Section || 
we express a number of mean-field quantities (such as 
the chemical potential /i, the gap parameter A, the 
pair-correlation length £ pa ir, and the bound-state en- 
ergy e ) in terms of the complete elliptic integrals of the 
first and second kinds with argument k 2 = (yl + ac* + 
x )/(2y / l + x 2 ) [where x Q = fi/A] ranging from (BE 
limit) to 1 (BCS limit). In Section III our analysis is ex- 
tended to quantities such as the phase coherence length 
Cp/iase and the sound velocity s, by considering the Gaus- 
sian fluctuations about the mean field. To keep the pre- 
sentation compact, the properties of the elliptic integrals 
used in our treatment are summarized in Appendix |a|. 
In Appendix [b| we report for the sake of comparison the 
solution for the two-dimensional case given previously by 
Ref. [ll]. 



II. ANALYTIC RESULTS AT THE MEAN-FIED 
LEVEL 

We consider the Hamiltonian (fi = 1) 

H = J2 Jdv^iv) Mr) 

+g /"drVjWVjWViWV^r) (2.1) 

where Vv (r) is the fermionic field operator with spin pro- 
jection a, m is the fermionic (effective) mass, and g — Vil 
is the strength of the short-range (contact) potential be- 
tween fermions with V < (Q being the volume occupied 
by the system). 
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At the zero temperature, the mean-field equations for 
the gap parameter A and the chemical potential fj, are ob- 
tained by a suitable decoupling of the Hamiltonian (2.1) 
and are given bycij 
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(2.2) 



(2.3) 



where k is the wave vector, N the total number of 
fermions, and 



A 2 , vl = 



1 



1 - 



in 

E k 



(2.4) 



Owing to our choice of a contact potential, Eq. < \2. 
diverges in the ultraviolet (both in two and three dimen- 
sions) and requires a suitable regularization. In three 
dimensions it is common practice to introduce-the scat- 
tering amplitude a s defined via the equationtltj 
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(2.5) 



where the divergent sum on the right-hand side of 
Eq.( |2.5D results in a finite value of a s b y letting V — > 
in a suitable way. Subtracting Eq. (2.5) from Eq. ( |2.2| ) 
we obtain 
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(2.6) 



which has to be solved simultaneously with Eq. (2.3) to 
determine A and fi as functions of a s . 

In three dimensions it is convenient to introduce the 
following dimensionless quantities 



*1± 
2m A 
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(2.7) 



with fc 
form 



|k|, and express Eqs. (|2.3|) and (|2.6| ) in the 



- - = -(2mA) 1 ' 2 I 1 {x ) 
a s 7T 

n= -i 7 (2TOA) 3/2 / 2 (a; 



(2.8a) 
(2.8b) 
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EZ ~ x 2 



c 00 

I\{x ) = i dxz 
Jo 

h(x ) = J dxx 2 ^- 



(2.9) 
(2.10) 



The integrals (2.9) and (2.10) were originally considered 
by Eaglesp who evaluated them numerically as functions 
of the crossover parameter x ranging from — oo (BE 
limit) to +00 (BCS limit ). We shall soon show that the 
integrals ( |2.9| ) and ( p. 10 ) can actually be evaluated in a 
closed form for all values of x . Note that l2(x ) > 
while Ii(x ) can take both signs. 

Before proceeding further, it is convenient to render 
Eqs. ( |2.8| ) dimensionless by introducing the Fermi en- 
ergy cf = k F /2m = (37r 2 n) 2 / 3 /2m where Uf is the Fermi 



momentum. In this way Eq. (2.8b) becomes 



A 



2/3 



and Eq. (2.8a) reduces to 



k F a s 
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h(x ) 



(2.11) 



(2.12) 



Note that the right-hand sides of E qs. ( |2.1l| ) and fl2.12j ) 
depend on x only. Equation (2.12) can thus be inverted 
to obtain function of kFd s ; from Eq. ( |2.1l[ ) and 

from n/tp = x Q A/eF one can then obtain the two pa- 
rameters A/ep and n/ep as functions of kpa a . 

Alternatively, one can use kp£, pa ir as the independent 
variable in the place of kpa s , where £, pa ir is the charac- 
teristic length for pair correlation given by the following 
expression in three dimensions at the mean-field level:Lj 
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where 



h{x ) = J dx^ 

f°° x 2 
h(x ) = dx— 



(2.14) 
(2.15) 



are two ad diti onal integrals expressed in terms of the 
quantities (2.7). Contrary to the integrals I\ and the 
integrals I3 and £4 are elementary and can be evaluated 
via the residues technique. One obtains: 



7T Xl(l + x\) 



h(x ) = -zXx 



(2.16) 
(2.17) 



where 



with the notation 
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T 2 



y/1 + X 2 + X a 



(2.18) 



Making use ofEq. ( PH ) we obtain eventually 
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(2.19) 



In this way Eq. (2.12) can be dropped in favor of Eq. 
(2.19), which can be inverted to obtain 

Xo as a function 

of kp^pair- 

Th ere re mains to evaluate the integr als I\ {x ) en terin g 
Eq. fl2.12j ) and h(x a ) entering Eqs. fl2.1l|) and ( ^.19[ ). 
To this end, we introduce the auxiliary integrals 



f°° x 2 
h(x ) = dx-j^ 

X £x 



h(x ) 
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El 



such that 



h{x ) = 2 (x I 6 (x ) - h(xo)) 
2 

h{x ) = - (x I 5 (x ) + h{x )) 



(2.20) 
(2.21) 

(2.22) 
(2.23) 



after integration by parts. The auxiliary integrals /s(x ) 
and Ie(x ) can, in turn, be expressed as linear com- 
binations of the complete elliptic integrals of the first 
[F(^,k)] and second [-E^k)] kinds. Reduction of a 
generic integral of the elliptic kind to the normal Legetu 
dre's form has been treated at length in the literature.Efl 
Here we proceed a s foll ows. Integration by parts gives 
for h{x ) [cf. Eq. (^2l|) 1: 



h{x ) 
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(2.24) 



where in the last line use has been made of the results of 
Appendix A and where (0 < k 2 < 1) 



K 2 = 



(1+ £2)1/2 



(2.25) 



with xi g iven by Eq. ( 2.18 ). For Iz{x ) we obtain instead 
[cf. Eq. (HI)]: 



h(x ) 



dx 



(a; 4 - 2x x 2 + x 2 a + l) 3 / 2 

1 . 7T . 



(2.26) 



where again use has been made of the results of Appendix 
A. 

In conc lusio n, we o btain f or th e quantities of interest 
[cf. Eqs. ( |2~Tll ), (flf ), and (p7lQ|) T : 



A 1 



£F (x I 5 (x ) + I 6 (x )) 2 / 3 
H [i A x 



e F Ae F (x I 5 (x ) + h{x )) 2/3 



1 



kpa s 

k F^pair 



4 (x I 6 (x ) - h(x )) 
7T (x I 5 (x ) + 7 6 (a; )) 1 /3 



l + xi 



4\ 1/2 



(x / 5 (x )+J 6 (a: )) 1 / 3 
(l + x 2 )^ 



(2.27) 
(2.28) 
(2.29) 

(2.30) 



with I$(x ) and h(x ) given by Eqs. ( |2.26 ) and ( 2.24 ), 
respectively. It is sometimes convenient to normalize 
the chemical potential, when negative, with respect to 
the bound-state energy e of the associated two-fermion 
problem. In the three-dimensional case, e G can 
pressed in terms of a s whenever a s > 0. One findi 



1 



mat 



so that [cf. Eq. fl2.29| )] 

e 2 32 (x I 6 (x ) -h{x )f 



(k F a s ) 2 7T 2 (x I 5 (x ) + I 6 (x )) 2/3 



(2.31) 



(2.32) 
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FIG . 1. A/eF vs kF^pair, obtained from Eqs. ( 2.27 ) and 
( j2.30|) . Full curve: exact solution; dashed curve: BCS ap- 
proximation, obtainedby including only the first two terms 
in Eqs. (A. 7) and (A.£); dotted curve: BE approxima tion , 
obtai ned by including only the first two terms in Eqs. (A.E) 
and (IO). 
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FIG. 2. h/ef for u > a nd u/ Un/2) f or u < us kF^pair, 
obtained from Eqs. ( 2.28 ), ( 2.30 ), and ( |2 . 3 2| ) . Conventions 
are as in Fig. [j]. 
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obtai ned from Eqs. ( 2.3C ) and 
(2.2£) is positive. Conventions 



e /tF VS ttF^pair 

(2.32) when a s given by Eq. 
are as in Fig. [j]. 



Numerical values of the complete elliptic integrals F 
and E have been extensively tabulated.^ 2 ! Otherwise, one 
may generate them with the required accuracy via Eqs. 
( [A.5| )-( pO| ). In this way, the desired values of h(x ) 
and Iq(x ) can be obtaine d for given x (or for given 
kptpair by inverting Eq. (pp|)). In Figs. [[[§ we re- 
port, respectively, the values of A/e F , fi/e F for /i > 
and /z/(e /2) for ji < 0, and e /e F vs k F ^ pa ir obtained 
by this procedure. Within numerical accuracy, all values 



coincide with those calculated by solving numeri cally the 
gap equation and the normalization condition (2.3) for 
the limiting case of a contact potential.t£l 

Besides the exact result (full curve), Figs. [l|j3] show 
for comparison two additional curves obtained by ap- 
proximating, respectively, the ell iptic integ rals F and E 
by the first two terms of Eqs. (A. 5) an d (|A.6| ) (d otted 
curve) and by the first two terms of Eqs. ( |A.7D and ( A.S) 
(dashed curve). In principle, these approximate results 
are expected to be reliable in the BE and BCS limits, 
in the order. Note, however, that the BE approximate 
result is surprisingly accurate on the BCS side of the 
crossover. 

In the BCS and BE limits the values of I§{x ) and 
Ie(x ) can be obtained by ret aini n g on ly a few signifi- 
cant terms in the expansions (A.5)-(A.8). In particular, 
in the BCS limit x a > 1 so that k 2 ~ 1 - \/(Ax 2 ) and 



(2.33) 



h(x ) ^ y^o 

h(x ) ~ \nx /(2^Q 



Then: 
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In the BE limit, on the other hand, x Q < and \x \ ^S> 1 
so that k 2 = 1/(4x1) and 



h{x ) 
h{x ) 



tt/(16|x | 3 / 2 ) 
7 r/(4|a ?0 | 1 / 2 ) 



(2.35) 



Then: 



A/e F 
fj./e F 
l/(k F a s ) 



[16/(37r)] 2 / 3 |ar | 1 /3 

-[16/(37r)} 2 / 3 \x \^ 
[16/(3tt)] 1 /3| ;Co |2/3 

2[16/(3tt)] 2 /3| Xo |4/3 
^[16/(3tt)] 



(2.36) 



-1/31 



-2/3 



The limiting BCS fl2.34| ) and BE ( |2.36D values coincide 
with those calculated previously by different methods Jl3 

We mention, finally, that another quantity which can 
be evaluated analytically at the mean-field level for all 
values of x a is the single-particle density of states. 

All the above results hold for the three-dimensional 
system. Analogous results for the two-dimensional sys- 
tem are straightforwardly expressed in terms of elemen- 
tary integrals and are reported for comparison in Ap- 
pendix pi 



III. ANALYTIC RESULTS AT THE GAUSSIAN 
LEVEL 

Besides the quantities of Section || defined at the 
mean-field level, additional quantities whose definition 
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requires the introduction of Gaussian fluctuations can 
also be expressed analytically in thr ee d imensions at zero 
temperature for the Hamiltonian (2.1), in terms of the 
complete elliptic integrals F and E for all values of the 
parameter x (i.e., following the evolution from BCS to 
BE). In particular, we shall consider the phase coherence 
length t; p hase (associated with the spatial fluctuations of 
the superconducting order parameter) and the sound ve- 



locity s (associated with the Goldstone mode of the bro- 
ken symmetry). 

The -matrix of the Gaussian fluctuations has 
elementM[L2 



r(q ' W) i%c) A(-q,-uj) 



(3.1) 



where uj is the frequency and (for a real order parameter) 



A(q,w) 



ly(j_ 
O ^ I 2E k 



2 2 



+ £k-q ~ U> - IT] E k + -Ek-q 



+ lu + irj 



(3.2) 



B(q,Uj) = ^ $>M-q«k-q [ E +E j 

k 



k-q -w — ir] E k + ,E k _ q + UJ + if] 



(3.3) 



at zero temperature. In these expressions, u 2 , — 1 — v 2 . 
and r] is a positive infinitesimal. To extract £ p hase and s 
we need consider the expansion of A(q, uj) and B(q,uj) 
for small values of |q| and w, such that 



2m 



where 



2A 



2v/A 2 + /i 2 



and uj <C uj* 



fi>0 
(i< 



The result is: 



A(q, uj) — a + a\ uj + a 2 q 2 + a 3 uj 2 
B(q, uj) = b a + b 2 q 2 + b 3 uj 2 + . .. 



with 



~ n ^ aeI 

k K 
_ 1 V- & 
k K 
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^ 32m 
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(3.4) 



(3.5) 



(3.6) 
(3-7) 



(3.8a) 



(3.8b) 



1 fCk(2£k - A 2 ) k 2 A 2 (8^ + 3A 2 



dm 



l + A 2 



(3.8c) 
(3.8d) 



and 



b = a 

52 = O 12 *k 



(3.9a) 



Q 32m 
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3^ A 2 k 2 A 2 (2£ 2 -3A 2 ) 
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fi ^ i6££ 

k K 



(3.9b) 
(3.9c) 



(d being the dimensionality). In particular, to determine 
iphase one has to consider the quantityE2l 



A(q, uj = 0) + B(q, w = 0) = 2a + (a 2 + & 2 )q 2 + 



from which 



P 2 

>>pha 



phase 



a 2 + b 2 
2a 



(3.10) 



(3.11) 



To determine the sound velocity s one has instead to 
consider the full determinant 



A(q,uj)A(q, -uj) ~ B{q,uj) 2 = 

2a (a 2 - b 2 )q 2 + [2a (a 3 - 63) - a\]uj 2 - 

which vanishes for uj = uj(q) = s\q\, with 

?2 _ 2a (q 2 - b 2 ) 
2a Q (b 3 — a 3 ) + a\ 



(3.12) 



(3.13) 



We are le ft wi th evaluating the integrals entering Eqs. 
( |3.11 ) and (3.13). In three dimensions we obtain: 
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(3.15) 
(3.16) 
(3.17) 
(3.18) 



where use h as be en ma de of the notation (2/7) and of 
the integrals Q2.20) and ( |2.2l|). T he four new integrals 
appearing in Eqs. ( 3.16 ) and fl3.17 ) can also be expressed 
as linear combinations of Is(x ) and lg(x Q ). Integrations 
by parts and simple manipulations lead to: 
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x 2 ix , 1 r, x A 



— / dx—^ 



(3.19) 
(3.20) 

(3.21) 
(3.22) 



In conclusion we obtain: 



a 2 + b 2 = 



1 1 /2m 



x i 2J 6 (x ) + 



(2tt) 2 12 V A 
(l + 4a;2) 



1/2 



3(1 + ^) 



[h(x )+x h{x )]\ (3.23) 



and 

a 2 - b 2 



1 1 /2m 



(2tt)2 3 V A 



1/2 



{h(x ) +x I 5 (x )} 



(3.24) 



Equation ( 3.11 ) thus reduces to 
e F 1 



A 12I 5 (x ) 



(1 + Ax 2 ) 

x { 2J 6 (a; ) + \J i °M l e (x ) + x I 5 (x )} 



3(l + xl)^°' ' "^"v-oyjj (3-25) 
with c_f / A given byEq. d2~27|), while Eq. (flf) becomes 
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1 A J 5 (a; )(J 6 (x ) + x / 5 (x )) 
3e F h{x ) 2 + h{x ) 2 



(3.26) 



where vf = kp/m. 

The expressions ( 3.25 ) and (|3.26 ) provide the desired 
analytic expressions of kp ^ v has e and s/^f for aZZ values 
of x . Using again Eqs. ( 2.24 ) and ( 2.26| ) and the tab- 
ulated values of the elliptic integrals, we report in Figs. 
H and|| the values of kp^ p h a se and s/vp, respectively, vs 

air • 



10 - 



1 - 




FI G. 4. kriphase vs kF^pair, obtained from Eqs. (2.3C) 
and (3.25). Conventions are as in Fig. jjj. 



To within numerical accur acy , we reproduce in this way 
the numerical results of Ref. [l(]for kp^ p h as e- Note, again, 
that the BE approximation (dotted curve) is surprisingly 
accurate even on the BCS side. 

In the BCS and BE limits one can use the approximate 
values ( 2.33] ) and (2.35), in the order, for the integrals 
l5(a; ) and Iq{x ). This gives: 



kp^phase 

s/v F 



x G /3 
1/V3 



(3.27) 
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in the BCS limit, and 



8/v F ~ (12ir\x \)-V 3 



(3.28) 



in the BE limit. Note that in the BE limit the product 
1 / * 



s£, P h a -s 



m \vf 



a $ 



l 

Am 



(3.29) 



is a constant and coincides with the Bogoliubov result 
(2ms) -1 for composite bosons with mass ms — 2m, thus 
confirming the general results established in Ref. [l(J for 
the mapping onto a bosonic system in the strong-coupling 
limit. 




0.1 1 10 100 

K f P « lr 

FIG . 5. s/vf vs kp^pair, obtained from Eqs. ( ^.3Cl|) and 
( 3.26 ). Conventions are as in Fig. [jj. 



collective mode associated with the symmetry breaking 
cannot be smaller than the size of a Cooper pair. 

The above results hold in three dimensions. In two di- 
mensions both £, p hase an d s can be expressed in terms of 
elementary integrals, as shown in Appendix p[ 



IV. CONCLUDING REMARKS 

In this paper we have provided the analytic solution 
of the crossover problem from BCS to BE in the three- 
dimensional case for a system of fermions interacting via 
an attractive contact interaction in free space, at the 
mean-field level and with the inclusion of Gaussian fluc- 
tuations. Although the assumptions required to obtain 
our analytic solution might be oversimplified in appli- 
cations to realistic systems, it could be interesting yet 
to compare our analytic solution with numerical calcu- 
lations describing more realistic cases. In particular, 
besides adopting a more sensible momentum-dependent 
form of the interaction potential, the free one-particle dis- 
persion relation ougth to be replaced by the actual band 
structure of the medium. 

In addition, a detailed description of the crossover 
problem from BCS to BE would require one to intro- 
duce already at the mean-field level the coupling with the 
charge degrees of freedom, whose effects are expected to 
be especially important in the crossover region of inter- 
est, intermediate between these two limits.t3 

In spite of these limitations, and considering the fact 
that analytic results for the crossover problem from BCS 
to BE were thus far limited to the two-dimensional case 
or to the two limits, our analytic solution is useful as it 
enables one to describe the crossover region in a compact 
way with very limited effort. 



An additional quantity which can be evaluated ana- 
lytically at the Gaussian level is the coefficient 7 of the 
quartic term in the dispersion relation 



c(q) 2 



2 2 
s q 



7 



Am 



(3.30) 



obtained by expanding the determinant ( 3.12 ) to higher 
order. Its expression in rather involved and will not be 
reported here. It is nonetheless interesting to note that 
in the BCS limit 



7 



256 
135 



(3.31) 



is negative (and large), so that the dispersion relation 

1 



W (q)< 



2 2 



1 - -q 2 £ 2 



(3.32) 



holds for |q| smaller then a critical value q c oc l/£ po ir- 
This implies that in the BCS limit the wavelength of the 
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APPENDIX A: RELEVANT PROPERTIES OF 
ELLIPTIC INTEGRALS 

In this Appendix we briefly review the properties of 
elliptic integrals that are relevant to our treatment. 

The elliptic integrals ofriiie-iirst and second kinds with 
modulus k are defined byllHlO 



F(a,K) = / dip—= 



1 



9 • 2 
k z sin ip 



E(a, k) 



dp 



1 9-2 

1 — k z sin ip 



(Al) 
(A2) 
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with k 2 < 1. They satisfy the properties (n integer), and are said to be complete when a = n/2. 

7r The complete elliptic integrals admit the following scries 

F(nn, K ) = 2nF(-, K ) (A3) representations: 

E(mr, k) = 2nE{^,n) (A4) 
I 



^^f{'-a) 2 « 2 -(^)^--(^) 2 ^r-] 

1 (, 4 1 \ , 2 , l 2 • 3 / 4 2 1 \ /4 



£(_ K ) = 1 + - hi U' + V— In «/+••■ (A 

V 2 ' 7 2l«' 1-2/ 2 2 -41k'1-2 3- 4/ V 



where «/ = VI — k 2 is known as the complementary 
modulus. The representations (A. 5) and (A. 6) are to be 
preferred when k 2 < 1; when k z ~ 1 and n' 2 -C 1 the 
representations (A. 7) and (A. 8) are to be preferred in- 
stead. 



Equations (2.24) and (2.26) of the text are obtained b 



adap ting tabulated results and using the properties ( 
and (|A4|) . We obtainEa 



dx - = — F(n, k) = -F(-,k) 

o Vx 4 + 26 2 a; 2 + a 4 2a v ; a V ; 



(A.9) 



f/.r; 



aE{n, k) 



F(tt,k) 



Vx 4 + 26 2 a; 2 + a 4 2(a 4 - 6 4 ) 4a(a 2 - 6 2 ) 



1 



a ,7T 
^ 4_ fe 4 S( 2' K) ~ 2a(a 2 - b 2 ) 



F(-,«) (A.10) 



where a 2 > b 2 > — oo a 2 > 0, an d k 2 = (a 2 - 6 2 )/ (2a 2 ) 
Com paris on of Eqs. (|A.9[) and (|A.10|) with Eqs. fl2.24|) 



and ( 2.26 ) leads us to identify 

6 2 = -x„ , 



a 4 = 1 + x 2 , 



y/1 + Z 2 , + z Q 
271+^1 
(All) 



and results eventually into the right-hand sides of Eqs. 
( pT24| ) and ( p726[ ). 



APPENDIX B: TWO-DIMENSIONAL CASE 

The analytic solution for the two-dimensional case has 
been already provided in Ref. [H]. Its derivation is sim- 
pler than for the three-dimensional counterpart obtained 



r 



in this paper, since it can be expressed in terms of ele- 
mentary integrals. For the sake of comparison, we report 
in this Appendix the two-dimensional solution on equal 
footing of the three-dimensional solution discussed in the 
text. 

Quite generally, in two dimensions the bound-state en- 
ergy e Q exists for any value of the interaction strength g. 
For the contact potential we are considering, however, 
the bound-state equation 



1 



k 2 /2 



m 
2^ 



dy 



1 



2y + e a /A 



(B.l) 



[with y = k 2 /(2mA)) needs to be suitably regularized by 
introducing an ultraviolet cutoff. This cutoff can, in turn, 
be re moved from further cons ider ation by combining Eq. 
( B.l ) with the gap equation (2.2), namely, 



g 2fl ^ 



1 



= — I dz 

A 2 4?r 



(1 



z 2)l/2 



(B.2) 



[with z = y — x and x Q given by Eq. (2.7)] which requires 



an analogous regularization. P erfor ming the elementary 
integrations in Eqs. (B.l) and (B^) one obtains 



The normalization condition (|2.3|) gives further 



(B.3) 



mA 
~2n 



dz[l- 



mA / — 

= — (x + ^T 



Since in the normal state n = kp/(2ir) 
(|B.4[) reads 



A 

Cp 



(B.4) 
mep/ir, Eq. 

(B.5) 
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Mult iplying at this point both sides of Eqs. (B.3) and 
(HH) yields 



X o 



X Q 



(B.6) 



Finally, the pair-correlation le ngth can be obtained 
from its definition [cf. Eq. (2.13) for the three- 
dimensional easel: 



e 2 



dk{k^i/El) 



dk(k/E%) 



2 Qy(ye v /E 6 v ) 

nA Qy{l/ED 



(B.7) 



with the dimensionless quantities 

f y = fc 2 /(2mA) , x = n/A 



Cy = V 



Ey - 



*2 



(B. 



[in the place of (2.7) for the three-dimensional case]. The 
integrals in Eq. (B.7) are again elementary and give: 



(hi: \ 2 — 1 eF 
{KFt,pair) — r 



— + arctan x 



(B.9) 



It is then clear that A/e F , n/e F = x A/e F , to/tF, and 
kF^pair can be expressed as functions of x a . [Note that no 
reference to the scattering amplitude a s has been given 
in two dimensions!. 



Alternatively, Eq. ( p.9\ j with e F / A given by Eq. ( |B.5| ) 
can be inverted to express x (as well as all other quan- 
tities) as a function of k F ^ pa i r . 

At t he Gaussi an l evel, all integrals entering the defini- 
tions (3.8) and (3.9) of the coefficients of the expansions 
(3.6) and (3.7) are elementary in two dimensions. We 
obtain: 



0,2 



a 
■b 2 



m yT+xf + x 
Btt y/TTxl 
m 1 



8tt A^l + x 2 
1 1 
48^A 



xf, + 3x 2 n + 1 



h - a s 



' (1+^)3/2 
1 tJ1+ x 2 Q + x 



16tt A 
m 1 yl + x 2 Q + x a 



i6tt a 2 yr 

These results give 



(B.10) 
(B.ll) 
(B.12) 

(B.13) 
(B.14) 



{kF^plu 



and 



6Ayi 



3xi 



1 



(1+^2)3/2 



(B.15) 



1 A 



(B.16) 



Note that, owing to Eq. (B^), (s/v F ) 2 = 1/2 is indepen- 
dent from x . On the other hand, for k F £, p i ia se we obtain 
in the two limits: 



k F ^phas 



X /V6 

1/V8 



BCS limit 
BE limit 



(B.17) 



thus confirming that the BE limit depends markedly 
on dimensionality and shows a peculiar behavior in two 
dimensions.t3 Note, however, that in the BE limit the 
product s^phase still coincides with the Bogoliubov re- 
sult (4m) _1 . 
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